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ABSTRACT 

Backscatter  of  acoustic  signals  generated  by  surfaces 
having  point  scattering  and  line  scattering  features  is 
studies  experimentally  in  air.   Signals  are  analyzed  for 
propagation  over  three  models:   a  smooth  plane  wedge,  a 
wedge  whose  surface  is  slightly  roughened  with  uniformly- 
packed  lead  shot,  and  one  randomly  roughened  with  gravel. 
Scattering  from  the  rough  surfaces  at  normal  incidence  is 
measured  to  estimate  the  statistical  parameters  that  describe 
the  surface  roughness.   The  returned  acoustic  energy  from 
the  composite  area  of  rough  surface  and  wedge  for  low  grazing 
angles  is  analyzed  to  determine  whether  the  identity  of  the 
two  individual  backscatter  effects  can  be  separated  in  their 
combined  spectrum.   The  near  grazing  backscattered  amplitudes 
are  much  larger  than  predicted  by  conventional  rough  surface 
scattering  theory  and  suggest  the  presence  of  a  backscattered 
boundary  wave. 
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I.   INTRODUCTION 

Diffraction  of  acoustic  energy  from  both  smooth  wedges 
and  rough  plane  surfaces  has  been  theoretically  investigated 
by  a  number  of  authors.   Many  at  sea  measurements  have  been 
reported  for  rough  surfaces.   Of  particular  interest  is  the 
theoretical  development  of  Biot  and  Tolstoy  [1]  for  impulse 
diffraction  of  point  source  radiation  by  plane  rigid  wedges, 
and  the  experimental  results  for  sea  surface  scattering, 
e.g.,  Clay  and  Medwin  [2]  and  Urick  [3]. 

Biot  and  Tolstoy  theory  has  been  tested  experimentally 
by  Bremhorst  [4]  and  Spaulding  [5]  for  simple  wedges  and 
complex  seamount  modeling,  respectively.   The  experimental 
results  showed  close  agreement  with  theory. 

This  study  continues  the  experimental  comparison  of 
diffraction  by  various  effects  in  the  specific  case  of 
backscatter.   The  approach  involved  determining  the  spectral 
relationship  for  the  smooth  wedge  and  a  rigid,  random  rough 
surface,  separately.   A  wedge  covered  with  a  rough  surface 
of  gravel  was  used  to  attempt  to  differentiate  between  the 
returned  signals  from  both  the  wedge  and  rough  surface 
effects.   Also,  a  wedge  covered  with  lead  shot  was  used  to 
present  a  uniformly-packed,  slightly  rough  surface.   The 
results  were  intended  to  show  how  the  combination  of  wedge 
slope  angle  and  surface  roughness  interact.   This  relationship 
would  ultimately  identify  a  range  of  these  characteristics 


where  one  signal  is  not  extremely  dominant  over  the  other, 
but  rather  the  two  signals  could  be  deconvolved  to  show  their 
separate  effects;  and  thereby  identifying  the  character  of 
the  source  of  backscatter. 


II.       RESEARCH    FACILITIES 

A.  ANECHOIC  CHAMBER 

The  experiments  were  conducted  in  the  Naval  Postgraduate 
School  anechoic  chamber  (air) ,  which  provides  excellent  iso- 
lation from  external  noise  sources,  absorption  of  internal 
sounds,  and  atmospheric  stability.   The  internal  design  of 
Fiberglas  wedges  allows  for  absorption  of  approximately  99% 
of  the  incident  sound  of  frequencies  greater  than  100  Hz. 

B.  DATA  PROCESSING  EQUIPMENT 

Data  acquisition  and  processing  were  accomplished  by 
using  a  digital  computer  system  composed  of  four  major  com- 
ponents interfaced  to  provide  high  speed  analog  to  digital 
conversion,  digital  processing,  and  data  printout.   The 
design  was  developed  by  the  Special  Projects  section  of 
Naval  Air  Development  Center  in  conjunction  with  Pinkerton 
Computer  Consultants,  Inc.,  of  Warminster,  Pennsylvania. 
The  major  components  are: 

1.  Interdata  Model  70  Computer 

This  minicomputer  is  a  digital  design  that  is  FORTRAN 
and  BASIC  programmable  with  a  6  4  thousand  byte  memory.   Addi- 
tionally, data  stored  on  floppy  discs  or  digital  tape 
cassettes  can  be  input  for  processing. 

2 .  Phoenix  Analog  to  Digital  Converters,  Model  ADC  912 
Two  ADC  912  analog  to  digital  converters  may  be  used 

separately  or  simultaneously.   Each  converter  is  a  high 


speed,  high  accuracy  device  capable  of  encoding  ±10  volt 
input  signals  into  twelve  binary  bits  of  data  with  a  reso- 
lution limit  of  one  part  in  4,095  at  a  maximum  rate  of  2 
microseconds  per  conversion.   It  measures  the  input  voltage 
relative  to  an  internal  precision  reference  voltage  with 
an  accuracy  of  ±0.025%  of  full  range.   The  ADC  912' s  will 
accommodate  a  typical  commutating  through-put  rate  of 
476,190  channels  per  second,  including  settling  time.   The 
sampling  frequency  is  set  in  a  stable  frequency  oscillator 
and  sent  to  the  converters  via  a  sampling  circuit.   As  used 
conservatively  at  NPS ,  the  maximum  sampling  frequency  is 
320,000  Hz. 

3.  Qrbis  Diskette  Drive,  Model  76S 

The  Model  76S  diskette  drive  is  a  peripheral  device 
interfaced  with  the  Interdata  70  Minicomputer  in  order  to 
read  and  write  on  an  industry  standard  floppy  disc.   The  dual 
drive  unit  provides  3.2  million  bits  of  data  with  a  data  rate 
of  250,000  bits  per  second  encoding  for  each  drive  and  is 
capable  of  a  6  millisecond  access  time  track-to-track. 

4 .  Texas  Instruments  Silent  Electronic  Data  Terminal, 
Model  7  33 

The  TI  733  consists  of  a  keyboard  used  as  a  program- 
ming input/output  control  device,  a  printer,  and  transmit/ 
receive  mechanism  for  the  peripheral  disc  drive  unit.   The 
overall  system  facilitates  rapid,  accurate  processing  of  any 
desired  type  of  analog  electrical  signal  and  was  used  pri- 
marily for  frequency  domain  analyses  using  standard  FFT 
algorithms . 
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C.   STANDARD  EQUIPMENT  LIST 

The  following  scientific  equipment  was  utilized  to 

conduct  the  experimental  work: 

Interface  Technology  Timing  Simulator/Word  Generator, 

Model  RS-648 

Wavetek  Arbitrary  Waveform  Generator,  Model  175 

Instruments,  Inc.,  Power  Amplifier,  Model  L2 

E&L  Instruments  DD-1  Digi  Designer 

Tektronix  Type  551  Oscilloscope  with  two  dual-trace 

plug-in  units 

General  Radio  Type  1312  Coherent  Decade  Oscillator 
Hewlett-Packard  Electronic  Counter,  Model  5223L 
Fluke  Digital  Multimeter,  Model  8000A 
Hewlett-Packard  DC  Power  Supply,  Model  721A  (two) 
Princeton  Applied  Research  Pre-Amplif ier ,  Model  113  (two) 
Krohn-Hite  Frequency  Filter,  Model  3550 
Bruel  and  Kjaer  Type  4145  Condenser  Microphone 
Bruel  and  Kjaer  Type  413  3  Condenser  Microphone 
Bruel  and  Kjaer  Microphone  Preamplifier,  Model  2619 
Bruel  and  Kjaer  Microphone  Power  Supply,  Model  2804 
The  equipment  diagram  for  the  experiment  is  shown  in 

Figure  3. 
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Figure  2 .   Typical  Rack  Mounting 
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III.   THEORY 

Scattering  can  be  generally  explained  as  the  phenomenon 
which  occurs  whenever  a  wave  front  encounters  an  obstacle, 
thus  producing  a  new  source  by  the  principle  of  Huygens  wave- 
lets.  With  the  exception  of  forward  scatter,  we  can  separate 
the  sound  waves  scattered  by  an  object  from  the  direct  sound 
wave  by  using  a  pinging  technique  and  observing  the  scattered 
signal  when  there  is  no  signal  from  the  source  [2] .   This 
experimental  work  involves  the  study  of  two  theories  for 
backscatter — that  for  which  sound  impinges  on  a  change  of 
slope  such  as  for  a  wedge  and  that  for  sound  scatter  from  a 
rough  surface. 

A.   WEDGE 

To  examine  the  former  theory,  consider  the  model  of  an 
infinite  wedge  which  is  bounded  by  rigid  plates  at  a  =  0 , 
6=9   and  intersecting  along  the  apex,  chosen  as  the  z-axis. 
The  region  is  occupied  by  a  homogeneous  compressible  fluid 
of  sound  speed,  c,  and  density,  p. 

To  determine  the  normal  modes  of  the  system,  it  is 
convenient  to  solve  the  acoustic  wave  equation  in  cylindri- 
cal coordinates,  in  terms  of  the  displacement  potential  $: 


2  2  2 

d    <t>         1  3  4>    1  9  $    3  (j>   _   1  3  4>  ,,, 

7?    ?^    7  7?    7?  '  7  7? 
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The  harmonic  solutions  to  this  equation  are  of  the  form 


♦<r,8,«,t)   =  e±iu6H(1'2)(Kr)e±iYZe±iut  (2) 

V 


where  v,  «,  and  y  are  separation  constants  and  may  be  thought 
of  as  components  of  the  wave  propagation  vector,  k,  in  the 
azimuthal,  radial,  and  z  directions,  respectively,  and 


2     ? 
— 2  ~  <  (3) 


Boundary  conditions  can  be  established  such  that  the 
sides  of  the  wedge  are  assumed  to  be  perfectly  rigid,  re- 
quiring the  displacements  normal  to  the  sides  to  vanish, 
and  can  be  expressed  mathematically  as, 

||  =  o    at   e  =  o,  e  =  ew  .        (4) 

Writing  the  first  term  of  the  RHS  of  equation  (2)  as 


±lv8  n     , 

e      =   cos  v9  +  sin  v 


and  applying  the  boundary  conditions  yields  the  solution 
to  equation  (1) , 


„  TT(1,2)   ±i>z   ±iwt  ,Cx 

=   cos  v  9  H   '    e   f   e  (5) 

n    v 
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where 


_   n  ^  I  c 

vn  "  TT  *  (6 

w 


At  r  =  0,  the  imaginary  part  of  the  Hankel  function 

becomes  infinite,  which  is  not  physically  realizable. 

(1  2) 
Therefore,  only  the  real  parts  of  H   '   ,  which  are  Bessel 

functions  of  the  first  kind,  will  be  kept. 

Biot  and  Tolstoy  [1]  develop  this  theory  by  assuming  an 

explosive  point  source  at  rn,  en,  and  z  =  0,  which  leads  to 

solutions  to  equation  (5)  that  are  symmetric  in  the  z-axis. 

Applying  the  method  of  normal  coordinates,  they  arrive  at 

the  solution, 


AH  f*  r-> 

r—     =      -   - —    >    cos    v    8    cos    v    9n    I       ,      ct>z       (7) 
3t  6TT   ■  n  n   0      n  — 

W   n 


where 


I       (  J   Ur)J   (<rn)  J[<(c2t2  -  z2)1/2]Kdic      (8) 
On      n 

For  a  detailed  look  at  the  mathematics  of  the  solution  tech- 
nique, see  Refs.  1  and  6.   A  transform  given  in  Appendix  X 
of  Ref.  6  may  be  used  to  reduce  I   to  more  physically  com- 
prehensible terms. 
Now  let 


fco    ■    ?[(r-ro'2+z2ll/2  (9) 
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lr,   L   .  2  J   2,1/2 
t'   =   -[  (r  +rQ)   +  z  ] 


where  t»  is  the  least  possible  time  required  for  a  sound 

pulse  to  travel  directly  from  a  source  to  receiver  and  t1 

is  the  time  required  to  travel  from  source  to  wedge  apex  to 

receiver.   Equation  (8)  will  have  three  solutions:   t  <  t~, 

t«  <  t  <  t '  ,  t '  <t.   For  t  <  t- ,  I   =  0  as  would  be  expected 
0  —   —        —  On  e 

and  there  is  no  signal.   For  tn    t  <  t1,  the  solution  des- 
cribes the  acoustic  field  for  the  direct  energy  (not  of 
interest) ,  reflected  energy,  and  any  backscattered  energy 
prior  to  the  signal  from  the  wedge  apex  (for  the  rough  sur- 
face case  considered).   For  t1    t,  the  solution  is  the  one 
of  interest  for  the  theory  applied  to  backscatter  from  wedges 
For  this  case, 


,  -v  y 

I =^r-r- —  sin  v  tt  e   n  (10) 

n       tt  r  rQ  sinh  y      n 


where 


2.2   ,  2  _,_   2^2, 
c  t  -  (r  +  rQ  +  z  ) 

y   =   arc  cosh  9    (11) 

z  r  rQ 


Using  this  solution,  equation  (7)  becomes 


-v    y 

H     =     — ri-r: —  y    cos    v    9ncosv    9    sin   v    rte      n      ,       (12) 

3t  tt8tt  rr.sinh  y     L  n    0  n  n 

WO  n 

t'    <    t 
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Rewriting  the  trigonometric  functions  in  terms  of  their 
exponential  identities  and  then  collecting  the  conjugate 
pairs,  equation  (12)  is  rewritten  as 


ii  =  _£ l e"TY/9w  (13) 

3t     4tt  6T7  r  rn  smh  y  v  o; 

WO       2 


-Try/6  -27ry/e  J 

l-2e       cos  (tt/8  )  (tt±9  ±eQ)  +e 

when  the  ±  signs  indicate  that  the  term  in  brackets  is  actualy 
a  sum  of  four  terms  corresponding  to  all  possible  combinations 
of  these  signs. 

For  the  unit  impulse  source  used  in  Refs.  1  and  6,  the 
displacement  potential  i>    is  related  to  the  acoustic  pressure 
by 

p   -   -  p  14  (14) 

3t 

If,  on  the  other  hand,  one  assumes  a  source  that  is  a  delta 
function  in  time  as  well  as  space,  the  acoustic  pressure 
resulting  from  the  backscattered  wave  alone  is  given  by 

-    p    c    1 e  W 

p  4tt    9tt  r  rn    sinh  y 

WO 

sin(7r/aw)  (iT±e±e0) 

*     S  =^7^  "2^/6        ]         (15) 

l-2e  cos(tt/9w)  (tt±8±60)   +e 
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where 


2^2   ,2     2    2X 

c  t  -  (r  +r   +  z  ) 

y   =   arc  cosh  


W 
}0 


2  r  rQ 


and 

=  fluid  filled  region  above  the  wedge 

=  angle  between  the  wedge  side  and  the  source 

=  angle  between  the  wedge  side  and  the  receiver 

r-.   =  distance  from  source  to  wedge  crest 

r   =  distance  from  wedge  crest  to  receiver 

z   =   distance  along  wedge  crest  measured  from 

point  of  intersection  of  least  time  travel 
path 

See  figure  4  for  illustration  of  geometry. 

The  backscattering  process  can  be  described  as  follows. 
A  point  source  at  a  distance  r~  from  the  wedge  crest  trans- 
mits a  spherical  wavefront.   The  wavefront  first  impinges 
the  crest  along  the  least  time  travel  path  at  time,  t  =  rn/c. 
The  wavefront  then  continues  to  arrive  at  progressively  later 
t  along  increasing  distance  z.   This  in  effect  produces  a 
continuous  series  of  point  sources  along  the  crest  which  re- 
radiate  .   Thus  the  crest  of  the  wedge  acts  as  a  line  source 
with  time  shading  along  the  line. 

Since  the  frequency  domain  solution  to  the  impulse  wave 
scattering  problem  is  more  familiar  to  the  applied  acoustician, 
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Figure  4 .   Typical  wedge  geometry 
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equation  (15)  has  been  transformed  to  the  frequency  domain 
by  Medwin  [7] ,  and  the  computed  theoretical  values  have  been 
compared  to  experimental  data. 

B.   ROUGH  SURFACE  SCATTER 

The  statistical  approach  to  the  scattering  theory  des- 
cribes the  random  surface  by  the  distribution  of  its  sur- 
face heights  about  a  mean  level  and  by  a  correlation  function 
that  relates  the  heights  at  various  points  on  the  surface. 
The  distribution  of  surface  heights  is  assumed  to  be  Gaussian 
with  zero  mean  and  standard  deviation,  a,  which  is  identically 
the  rms  surface  height.   The  correlation  function  is  assumed 
to  be 


C(t)   =   exp  -(t/L)2  (16) 


where  L  is  the  correlation  distance  of  the  surface  and  t  is 
the  distance  between  any  two  points  whose  correlation  is 
desired. 

The  derivation  of  this  theory  is  given  by  Beckmann  and 
Spizzichino  [8]  for  electromagnetic  radiation.   However, 
acoustic  scattering  is  analagous  to  electromagnetic,  where 
the  pressure  field  is  developed  by  using  the  Helmholtz  inte- 
gral and  applying  boundary  conditions  in  conjunction  with 
the  Kirchhoff  approximation  [2] .   The  further  assumption 
that  the  incident  pressure  is  constant  over  a  rectangular 
insonified  area  and  zero  elsewhere  is  imposed  to  yield 
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a  formula  for  the  scattered  pressure  from  a  random  rough 
surface, 


P2   =   p20  I  /  /  elV'r  dxdy  (17) 


P20   =   specular  scatter  from  smooth  surface 

i  k  p. Q  A  cos  9, 
p20   =  fTZ ^expdkr)  (18) 


p. Q  =  incident  pressure 

9.  =  angle  of  incidence 

9_  =  angle  of  scatter  to  receiver 

r~  =  range  from  surface  to  receiver 


1  +  cos  9,  cos  9«  -  sin  e ..  sin  9_  cos  <j> 

cos  91 (cos  e,  + cos  9T)  v 


v    =   propagation  wave  vector,  a  function  of 
geometry. 


Since  p_  is  a  complex  quantity,  the  mean  square  scattered 
pressure  <p9p9*>  must  be  calculated  in  order  to  introduce 
the  scattered  intensity.   This  then  is  related  to  the  mean 
and  variance  of  the  scattered  pressure  by  the  expression 

<p  p  *>   =   <P0>*<P0*>  +  Var[p9]  .  (20) 
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The  mean  scattered  pressure  is  equivalent  to  the  coherent 
component  of  the  specularly  scattered  pressure  field  since 
it  has  constant  phase  with  respect  to  the  smooth  reflected 
pressure. 

<P2>   =   P90  F  exp(-g/2)  sine  v  X  sine  v  Y    (21) 


where 


sin  v  X 

sine  v  X  =   tz — 

x        v  X 
x 


and  g  is  defined  as  the  roughness  parameter. 


2  2  2 

g   =   k  o     (cos  9,  +  cos  9?)   .  (22) 


The  variance  is  derived  by  assuming  surface  roughness  to 
be  isotropic, 


2   2  2 

*  F   L   exp(-g)  p9 

Var[p2]   =   j 20_ 


00   m  -  (v   +  v   )  L 

.  I      2_^  exp[ - — j-2- ]      (23) 

S  mm!  ^        4m 
m=l 


and  this  equates  to  the  incoherent  component  of  the  scattered 
pressure  field.   A  more  detailed  analysis  of  the  mathematics 
in  this  section  can  be  found  in  Ref.  8. 
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For  the  specular  case,  the  scattered  intensity  reduces 
to 

<p  p  *>  2  m 

-~T     '      exP("^  +  ^[exp(-g)  I      2_]    (24) 
|P20I  m=l 

The  expression  in  brackets  is  defined  as  S(g),  and  approaches 
a  constant  linear  slope  at  very  low  or  very  high  values  of 
roughness , 

S  (g)   ->   g    for   g  <  .  1 
S(g)   ->   1/g    for   g  >  10  . 

For  normal  incidence,  the  relative  scattered  intensity 

shows  a  relationship  to  frequency  squared  through  e  "   for 

2    2  2 
g  <  1  and  approaches  an  asymptotic  value  of  AL  /2a  r   for 

g  >  10. 

In  the  backscatter  case,  a  surface  backscattering  factor 

is  defined  as 


<PoP0*>    r. 


2 


2^2       2 

where  p.  is  the  incident  pressure  at  the  surface.   For  a  rela- 
tively smooth  surface,  g  <<  1,  the  intensity  becomes 

2  2  2 
16  tt  L  a  p2Q2        2  2    2 

<p  p  *>   =   — = 5 exp(-k  L  sin  Q.   -  g)  .  (26) 

Z    Z  AZ  cos  9   A  X 
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and  the  backscattering  factor  is 


S 


3  2  2 
16  ^JL  a       ,  ,  2,2   .2 


B         X4 


exp(-k  L   sin  e,  -g)  .      (27) 


For  the  rough  surface,  g  >>  1,  the  coherent  term  is  considered 
negligible  so  that  the  variance  equals  the  intensity 

3    2  2 

TrF  p20  L  2     2   2 

<P2P2*>   =   VartP2]   =   Ag" exp[-(vx  +  vy  )  L  /4g] 


(28) 


29) 


and  the  backscattering  factor  becomes 


where 


SD   =   (8tte2  cos461)~1  exp[-tan   6  /2s  ] 
B  J.  J- 


e        =  -—y-  =   mean  square  slope  .  (30) 
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IV.   EXPERIMENTAL  DESIGN  AND  PROCEDURES 

A.   DESIGN  OF  THE  MODELS 

The  choice  of  models  to  experimentally  test  the  theory  was 
a  progression  from  a  simple  wedge  to  more  complex  roughened 
surfaces  and  finally  to  a  combination  of  these  two  features. 
It  was  necessary  that  the  models  appear  to  the  source  and 
receiver  as  being  spatially  infinite  compared  to  the  wave- 
length of  the  signals.   By  using  sufficiently  large  model 
dimensions  and  a  pulsed  signal  technique,  interference  from 
edge  diffraction  was  avoided  and  the  surface  could  be  con- 
sidered infinite  for  the  frequency  band^of  the  experiment. 
Additionally,  the  models  were  required  to  approximate  a  per- 
fectly rigid  condition.   Aluminum  was  selected  for  the  build- 
ing material  due  to  its  availability,  machinability ,  rigidity, 
and  reflection  coefficient  in  air  of  almost  unity. 

1.   Smooth  Wedge 

The  wedge  fabricated  for  Spaulding's  [5]  experimental 
work  was  available  and  appropriate  for  this  part  of  the  study. 
The  dimensions — two  pieces  of  1/4"  thick  aluminum,  each  1.52 
m.  (5  ft.)  long  by  0.6  m.  (2  ft.)  wide,  joined  at  the  crest 
forming  an  interior  angle  of  152° — were  dictated  by  the  size 
of  an  aluminum  section  available,  which  satisfied  the  infinite 
surface  criterion.   The  symmetrical  angle  of  14°  with  the 
horizontal  was  chosen  to  approximate  the  average  upslope 
and  downs lope  of  the  seamount  modeled  by  Spaulding.   For 
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this  study  of  backscatter,  the  wedge  presented  a  depression 
angle  of  28°  as  the  sound  wave  traveled  up  the  slope  and 
encountered  the  effect  from  the  crest. 
2 .   Rough  Surface 

The  rough  wedge  design  required  a  random  surface 
whose  rms  height  and  rms  slope  of  roughness  were  known 
parameters  for  comparison  to  theory.   Since  the  final  model 
was  to  be  a  combination  of  a  wedge  and  rough  surface,  the 
rough  surface  was  prepared  on  a  model  of  identical  dimensions 
as  those  for  the  smooth  wedge.   Due  to  weight  and  cost, 
aluminum  of  3/16"  thickness  was  chosen,  yet  the  massiveness 
was  maintained  by  the  addition  of  lead  shot  as  the  roughening 
element.   Since  the  wedge  had  a  dual  use  with  another  experi- 
mental study,  lead  shot  was  chosen  to  be  compatible  for  both 
experiments  and  was  distributed  uniformly  in  a  dense  packing 
pattern  over  the  surface.   It  was  intended  that  the  frequency 
range  of  the  study  be  kh  <  1  where  h  is  the  separation  of 
the  shot  so  that  the  regularity  of  the  uniform  packing  would 
not  affect  the  backscatter. 

However,  for  the  wavelengths  available  (due  to  the 
kh  <  1  limitation) ,  the  size  of  the  lead  shot  and  the  uniform 
packing  presented  a  surface  that  appeared  "acoustically 
smooth"  i.e.,  g  <  1.   Furthermore,  since  the  theory  for  a 
random  rough  surface  assumed  Gaussian  distribution  for  the 
roughness,  a  different  roughening  element  was  required  to 
approximate  this  distribution.   The  wedge-shaped  base  remained 
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the  same  with  gravel  applied  to  roughen  the  surface.   The 
statistical  quantities  of  this  surface  were  determined  from 
height  measurements  by  Lt.  E.  Jordan. 

B.   SOURCE/RECEIVER  SELECTION 

Bremhorst  [4]  conducted  extensive  experimental  testing 
to  determine  the  most  suitable  source  and  receiver  to  satisfy 
the  requirements  of  size,  frequency  response,  acoustic  out- 
put, and  receiver  sensitivity.   After  a  compromise  of  these 
parameters  for  the  best  experimental  combination,  Bruel  and 
Kjaer  (B&K)  condenser  microphones  were  chosen  as  both  source 
and  receiver. 

To  examine  the  Biot-Tolstoy  theory,  it  is  desirable  to 
approximate  a  point  source  to  transmit  spherically  symmetri- 
cal waves.   This  condition  requires  that  the  source  dimen- 
sions be  smaller  than  the  wavelength  of  the  transmitted 
signal  (ka  <<  1,  where  k  is  the  wave  number  and  a  is  the 
radius  of  the  source) .   To  approach  the  assumption  of  a  point 
source,  the  smallest  possible  microphone  should  be  used  as 
a  source;  however,  the  limiting  factor  for  size  is  the  need 
for  a  usable  acoustic  pressure  output  over  the  frequency 
range  of  interest.   This  was  especially  critical  when  trans- 
mitting over  the  rough  surface,  as  the  effects  of  secondary 
scatter  added  to  the  transmission  loss.   Although  the  B&K 
condenser  microphones  operate  as  piston  radiators,  they 
follow  the  fundamental  property  of  radiating  a  spherical 
diverging  wave  but  with  a  directivity  factor  included.   The 
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directivity  term  could  be  ignored  as  long  as  the  sound 
field  was  essentially  constant  over  the  insonified  area  of 
interest,  i.e.,  on  or  near  the  acoustic  axis  of  the 
microphone. 

The  receiver  selection  involved  the  trade-off  between 
the  smallest  possible  microphone  for  good  high  frequency 
response  and  the  largest  receiver  to  maintain  a  sensitivity 
level  sufficient  to  detect  the  weak,  returned  signals.   The 
optimum  combination  was  using  the  1"  B&K  4145  as  the  source 
and  the  1/2"  B&K  4133  as  the  receiver. 

Another  important  consideration  was  that  the  pulse  length 
of  the  single-ping  signal  determined  the  resolution  of  the 
surface  roughness  that  could  be  recognized.   The  B&K  conden- 
ser microphone  exhibited  the  disadvantage  of  transforming 
a  single-ping  waveform  into  a  repetitive  periodic  waveform 
decaying  with  time,  or  a  "ringing"  effect.   The  result  was 
a  larger  area  of  insonif ication  than  desired  because  of  the 
increased  time  of  the  output  signal. 

This  study  embarked  on  the  task  of  constructing  and  test- 
ing an  alternate  microphone  that  could  generate  the  required 
acoustic  power  while  maintaining  receiver  sensitivity,  yet 
reduce  or  remove  any  ringing  effect.   A  solid  dielectric 
Mylar  transducer  was  designed  using  the  guidelines  of  Ref.  9. 
Four  transducers  were  fabricated — two  3"  and  two  1",  with 
each  size  having  different  degrees  of  roughness  of  the  back- 
plate.   It  was  projected  that  if  these  designs  decreased  the 
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ringing  effect,  a  newer  material  of  polyvinylidene  fluoride 
(PVF„)  could  replace  the  Mylar  surface  to  further  enhance 
the  output  capability  and  the  receiver  sensitivity.   This 
material  has  the  advantage  of  being  prepolarized  to  allow 
more  electrical  signal  to  either  generate  or  detect  the  acous- 
tic signal. 

Testing  indicated  that  the  Mylar  and  PVF-  dielectric 
transducers  also  showed  a  ringing  effect,  so  that  three 
additional  tests  were  conducted  to  compare  the  amplitudes 
and  duration  of  received  signal  of  all  the  transducers.   A 
common  source  was  used  to  radiate  to  the  various  receivers, 
and  likewise  a  common  receiver  was  used  to  compare  the 
transducers  as  sources.   The  sources  and  receivers  were 
facing  each  other  on  their  respective  acoustic  axes.   The 
third  test  used  a  common  receiver  with  the  various  sources, 
but  compared  the  reflected  signal  due  to  normal  incidence 
on  a  rigid  plate  (Table  1) . 

A3"  Mylar  transducer  was  used  as  a  common  source  with 
the  resulting  differences  in  receiver  characteristics  being 
greater  than  any  of  the  other  tests.   The  1"  B&K  microphone 
(source)  had  a  signal  strength  23-29  dB  greater  than  any  of 
the  dielectric  transducers,  but  demonstrated  a  ringing  of 
600  microseconds  that  was  from  two  to  four  times  longer  than 
the  others . 

With  an  1/8"  B&K  microphone  as  a  common  receiver,  the 
sources  were  more  similar  in  their  results.   The  3"  Mylar 
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Table  1 

Transducer  Comparison — All  values  determined 
from  oscilloscope  photos 


A.   3"  Mylar  source 

RCVR  AMPLITUDE  [dB  re  lv] 

1/8"  B&K 
1/2"  B&K 
1"  B&K 
1"  PVF 
1"  Mylar 
3"  Mylar 


-18, 

.4 

-  1, 

.4 

0, 

.8 

-22, 

.0 

-28, 

.0 

-24, 

.4 

DURATION  [usee] 


150 

300 

600 

300 

200 

150 

B.   1/8"  B&K  receiver 

SOURCE  AMPLITUDE  [dB  re  lv] 


DURATION  [  y sec] 


1/2"  B&K 
1"  B&K 
1"  PVF2 
1"  Mylar 
3"  Mylar 


-16.5 
-14.0 
-16.5 
-14.0 
-10.5 


50 
100 
50 
40 
25 


C.   Reflection  off  smooth,  rigid  plate  at  normal  incidence; 
1/2"  B&K  receiver 

IB  re  lv]     DURATION  [ysec] 


180 
160 
120 
100 


SOURCE 

AMPLITUDE 

1" 

B&K 

-  3.1 

1" 

PVF2 

-18.8 

1" 

Mylar 

-18.4 

3" 

Mylar 

-26.0 

Note:   "Duration"  is  measured  from  onset  of  signal  to 
time  when  noise  is  greater  than  signal. 
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transducer  had  a  6  dB  stronger  signal  than  the  1/2"  B&K 
and  the  1"  PVF„  and  was  4  dB  stronger  than  the  1"  B&K  and  the 
1"  Mylar.   Only  a  small  variation  in  the  ringing  existed 
where  the  3"  Mylar  signal  duration  was  25  microseconds  as 
compared  to  the  1"  B&K  of  100  microseconds. 

The  reflected  signal  indicated  that  the  1"  B&K  source 
signal  exceeded  those  of  the  dielectrics  by  15  to  23  dB, 
yet  the  ringing  of  180  microseconds  duration  was  only  slightly 
more  than  the  best  case  of  100  microseconds  for  the  3"  Mylar. 

These  results  verified  the  known  fact  of  poor  efficiency 
for  dielectric  transducers.   Although  the  B&K  microphones 
were  limited  in  the  amount  of  applied  voltage  when  used  as 
a  source,  twice  as  much  voltage  applied  to  the  dielectric 
transducers  still  yielded  a  much  weaker  signal.   The  differ- 
ence of  the  ringing  characteristic  was  considered  less  signi- 
ficant as  a  variation  of  100  microseconds  produces  a  3.45  cm. 
spatial  resolution  change.   Also  an  inherent  discrepancy 
arose  during  the  testing.   The  electret  transducers  appeared 
to  produce  results  that  decayed  with  repeated  trials,  indi- 
cating that  the  stability  of  the  laboratory  fabrication  of 
the  transducers  was  not  on  the  same  scale  as  B&K  commercial 
production. 

The  conclusion  was  to  abandon  use  of  a  dielectric  trans- 
ducer and  revert  to  the  previous  findings--a  bistatic  com- 
bination of  a  1"  B&K  source  and  1/2"  B&K  receiver  [Ref .  10] . 
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C.   SIGNAL  PROCESSING 
1.   Source  Signal 

In  order  to  obtain  spectral  analysis  of  a  broad 
frequency  range,  the  transmitted  acoustic  signal  waveform 
must  be  composed  of  many  frequencies.   Ideally,  an  impulse 
would  provide  equal  energy  at  all  frequencies.   However, 
a  triangular  waveform  approximates  an  impulse  on  its  posi- 
tive going  pulse  and  provides  signals  at  all  frequencies 
with  harmonics  of  the  fundamental  predominant.   The  harmonic 
amplitudes  decay  at  a  rate  of  6  dB  per  octave. 

The  Wavetek  Model  175  Arbitrary  Waveform  Generator 
was  used  to  generate  the  signal  from  stored  waveforms  or  by 
programming  any  desired  shape.   A  rise  time  of  500  nano- 
seconds of  the  triangular  waveform  allowed  a  close  approxi- 
mation to  the  infinite  slope  of  an  impulse.   The  waveform 
generator  was  triggered  by  the  Interface  Technology  Model 
RS-648  timing  simulator  to  provide  a  single  pulse  mode.   This 
method  presented  the  ability  to  separate  the  deisred  scattered 
signal  from  other  interference  for  each  cycle  of  the  output 
waveform. 

The  pulsed  output  was  amplified  by  the  Instruments, 
Inc.,  Model  L-2  Power  Amplifier  to  the  maximum  allowable 
voltage  for  the  B&K  microphones.   For  the  1"  B&K,  a  total 
AC  + DC  voltage  of  250  V  is  the  limit.   The  best  combination 
of  polarizing  voltage  and  AC  driving  voltage  was  150  V  and 

i 

100  V  (peak  to  peak) ,  respectively. 
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2  .   Received  Signal 

The  received  acoustic  energy  at  the  microphone  was 
then  amplified,  bandpass  filtered  to  eliminate  low  frequency 
noise  and  prevent  aliasing,  and  amplified  again  to  approach 
the  maximum  input  limit  of  the  A/D  converter,  ±10  volts. 
The  amplifiers  used  were  PAR  Model  113  preamps  due  to  their 
very  low  self-noise  characteristics. 

To  separate  the  diffracted  signal  of  interest  from 
the  stronger,  unwanted  signals  of  direct  path  travel,  re- 
flected travel  from  the  surfaces,  and  diffraction  from  the 
edges;  a  sampling  circuit  was  used  to  select  a  window  coin- 
ciding at  the  time  of  the  diffracted  analog  signal.   The 
Interface  Technology  timing  simulator  accomplished  the  gating 
for  the  sampling  frequency  to  100  nanosecond  accuracy. 

The  sampling  circuit  was  built  with  IC  devices,  asso- 
ciated power  supplies,  and  an  oscillator.   A  General  Radio 
Model  1312  coherent  decade  oscillator  delivered  the  sampling 
frequency,  since  its  stability  characteristics  provided  the 
A/D  converter  with  a  constant  frequency  which  allowed  in- 
creased accuracy  and  repeatability  of  the  data  analysis.   A 
schematic  diagram  of  the  circuit  is  illustrated  in  Figure  8. 
The  circuit  operates  in  the  following  manner.   The  decade 
oscillator  porduces  a  constant  frequency  sine  wave,  which 
is  then  converted  to  the  square  wave  form  required  by  a  A/D 
converter  by  a  LM  710CN  voltage  comparator.   The  stable 
square  wave  output  of  the  voltage  comparator  is  combined 
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Figure  8 .   Sampling  circuit 
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with  the  timing  simulator  pulse  of  duration  for  the  desired 
aperture  in  a  SW  740 8N  AND  gate.   The  resulting  output  of 
the  AND  gate  is  a  constant  sampling  frequency  to  the  A/D 
converter  only  for  the  duration  of  the  window.   The  input 
pulses  to  the  A/D  converter  keyed  the  device  to  sample  the 
analog  signal  being  routed  to  the  computer  through  the  A/D 
converter. 

Within  a  sampling  aperture,  the  number  of  samples 
collected  by  the  A/D  converter  was  set  to  be  a  power  of  two; 
i.e.,  64,  128,  256,  etc.,  as  required  by  the  FOURONE  FFT 
algorithm,  which  performed  the  transform  for  the  frequency 
domain  analysis.   The  sampling  frequency  was  set  at  an  inte- 
gral multiple  of  the  frequency  of  the  transmitted  signal. 
This  relationship  could  be  controlled  by  the  decade  oscilla- 
tor and  the  Wavetek  175,  thereby  minimizing  any  truncation 
of  the  sampled  waveform  and  eliminating  sidelobes.   The  fre- 
quency resolution  was  determined  by  the  following  formula, 

sampling  frequency      -  .  . . 

fi-r- — s a t — *-  =   frequency  resolution 

number  of  samples         ^    a 

However,  the  limiting  value  of  frequency  resolution  was  fixed 
by  the  duration  of  the  sampling  window. 

The  experiment  was  conducted  with  a  5  kHz  transmitted 
signal  and  a  200  microsecond  sampling  window,  which  likewise 
set  a  5  kHz  frequency  resolution  in  the  output  data.   With 
this  resolution  determined-,  then  a  sampling  frequency  of  320 
kHz  would  require  64  samples  of  data  within  the  sampling 
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aperture.   Following  one  cycle  of  signal  transmission  and 
received  signal  processing,  a  20  millisecond  delay  was 
incorporated  in  the  timing  sequence  to  allow  for  complete 
decay  of  any  reverberation  before  the  next  ping  was  triggered 

To  improve  the  signal-to-noise  ratio,  an  averaging 
routine  within  the  computer  software  was  performed  on  the 
digitized  data.   An  ensemble  of  as  many  as  9,999  sampling 
windows  could  be  selected  for  averaging  prior  to  Fourier 
analysis.   The  signal-to-noise  ratio  could  be  increased  by 
10  log  /N,  where  N  is  the  number  of  averaging  runs.   Con- 
sidering the  trade-off  of  computer  run  time  versus  signal 
gain,  1,000  blocks  of  data  were  averaged  to  yield  a  15  dB 
improvement.   This  averaged  data  was  then  used  by  the  FFT 
algorithm,  so  that  the  computer  output  presented  intensity 
by  spectral  distribution. 
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V.   EXPERIMENTAL  RESULTS  AND  ANALYSIS 

A.   BEAM  PATTERN 

Prior  to  determining  the  reradiated  energy  from  the  wedge 
or  the  rough  surface,  a  beam  pattern  for  the  adjacent  source- 
receiver  combination  was  required  to  use  as  the  reference 
pressure  at  a  distance  of  one  meter.   The  source-receiver 
combination  (Fig.  9)  was  aligned  to  transmit  at  normal  inci- 
dence toward  a  rigid,  smooth  plate.   The  source-receiver 

o 

set-up  was  then  rotated  through  ±90  ,  so  that  the  received 
signal  described  the  beam  pattern  for  the  source-receiver 
combination.   This  procedure  was  performed  for  two  arrange- 
ments that  are  called  the  symmetrical  and  asymmetrical  cases. 
In  the  symmetrical  case,  the  source  and  receiver  are  aligned 
relative  to  each  other  vertically  (parallel)  with  respect 
to  the  axis  of  rotation,  and  in  the  asymmetrical  case,  they 
are  horizontal  (perpendicular)  to  the  axis.   As  the  combination 
is  rotated  in  azimuth,  both  the  source  and  receiver  are  the 
same  distance  from  the  rigid  reflector  for  the  symmetrical 
situation.   However,  in  the  asymmetrical  arrangement,  either 
the  source  or  receiver  is  closer  to  the  rigid  plate  than  the 
other  depending  upon  the  direction  of  rotation. 

The  resulting  beam  patterns  for  the  two  situations  are 
presented  in  Figures  10  through  14  for  various  frequencies. 

Interestingly,  the  two  patterns  are  quite  similar  regard- 
less of  the  relative  position  of  the  source  and  receiver  to 
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Figure  9 .   Beam  pattern  source-receiver  arrangement 
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Figure  10.   Beam  pattern  for  10  kHz 
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Figure  11.   Beam  pattern  for  30  kHz 
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Figure  12.   Beam  pattern  for  50  kHz 
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Figure  13.   Beam  pattern  for  70  kHz 
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Figure  14.   Beam  pattern  for  100  kHz 
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the  reflecting  plate.   This  can  be  accounted  for  by  the  fact 
that  the  200  microsecond  aperture  of  the  received  signal 
processing  is  long  enough  to  include  the  diffraction  results 
from  the  adjacent  microphone  in  either  spatial  arrangement. 
The  dissimilarity  between  the  right  and  left  halves  of  the 
pattern  is  likewise  consistent  throughout  the  frequency 
range  and  results  from  the  diffraction  effects  of  the  source- 
receiver  combination  and  the  mounting  set-up. 

B.   SMOOTH  WEDGE 

To  measure  the  backscatter  from  the  crest  of  the  wedge, 
the  source-receiver  acoustic  axis  was  aimed  at  the  crest 
(z-axis  of  wedge) .   The  position  of  the  source-receiver 
relative  to  the  wedge  was  varied  to  compare  to  the  theory  for 
the  respective  geometry.   The  position  was  defined  by  the 
slant  range  from  the  crest  to  the  source-receiver  and  by  the 
angle  between  the  acoustic  axis  and  the  wedge.   The  geometry 
of  the  experiment  was  limited  by  the  interference  of  unwanted 
signals  when  the  crest  backscatter  could  not  be  individually 
gated  from  the  other  returns.   Figure  15  shows  an  oscillo- 
scope trace  of  the  received  signal  with  the  three  returns 
separated  in  time.   The  initial,  large  amplitude  (clipped 
in  picture)  signal  is  the  direct  wave  which  is  received 
immediately.   The  reflected  signal  due  to  normal  incidence 
of  the  side  lobe  appears  at  approximately  820  microseconds, 
and  the  backscattered  return  from  the  crest  arrives  at 
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Figure  15.   Typical  oscilloscope  trace  of  direct, 
reflected,  and  diffracted  signals 
from  smooth  wedge 
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approximately  1.5  milliseconds.   The  square  pulse  above  the 
received  signal  shows  the  sampling  aperture  for  gating  the 
return  to  be  analyzed. 

The  geometry  of  Figure  4  and  the  picture  in  Figure  15 
show  how  the  source-receiver  position  limitations  are 
imposed.   If  the  source-receiver  is  too  close  to  the  crest, 
the  backscatter  will  be  masked  by  the  stronger  direct  and 
reflected  returns.   If  the  range  is  kept  constant,  but  the 
angle  is  increased,  the  reflected  signal  is  delayed  in  time 
and  overshadows  the  backscatter  again.   If  the  range  is  in- 
creased to  more  than  half  the  length  of  the  side  of  the  wedge, 
diffraction  from  the  back  edge  interferes  with  the  diffraction 
of  interest  from  the  crest.   To  determine  the  position  limi- 
tations numerically,  data  collection  was  conducted  at  angles 
of  0°,  10°,  20°,  30°  at  both  25  cm.  and  30  cm.   Figures  16 
through  23  show  the  normalized,  backscattered  intensity 
plotted  versus  frequency.   It  is  noted  that  the  backscatter 
for  a  wedge  can  be  defined  by  the  quantity  of  backscattered 
signal  strength  (BSS) , 

BSSW   =   10  log[-|£f-£] 
w  j.q  rQ  a 

where  I   is  the  intensity  at  one  meter.   However,  since  it 

is  the  intensity  ratio  which  is  calculated  in  the  signal 

processing,  this  is  the  quantity  used  in  the  comparison  with 

theory.   In  theory,  I    is  inversely  proportional  to  frequency 

Bo 

for  a  wedge  at  high  frequencies  (Ref.  7). 
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The  received  signals  were  time  averaged  1000  times  in  a 
200  microsecond  sampling  aperture  to  obtain  a  signal-to- 
noise  ratio  (S/N)  ranging  from  5  to  30  dB  throughout  the 
frequency  spectrum;  Figure  24  shows  the  relative  S/N  for  the 
three  surfaces.   The  noise  was  measured  with  both  source 
and  microphone  on  but  with  a  data  window  corresponding  to 
a  region  of  nonscatter  in  the  anechoic  chamber. 

The  results  of  Figs.  16-23  indicate  a  major  consistent 
deviation  from  theory  for  the  frequency  range  above  50  kHz. 
This  can  be  explained  by  reverting  to  the  beam  pattern  dia- 
grams.  At  the  ranges  of  the  experiment,  the  200  microsecond 
processing  window  for  the  returned  signal  subtends  a  half- 
angle  of  28°.   This  angle  of  insonified  area  should  be  within 
the  main  lobe  of  the  source-receiver  pattern  to  insure  symmetry 
and  approximate  a  point  source.   The  plot  for  50  kHz  approaches 
this  critical  condition  where  the  rapid  fall-off  from  the 
main  lobe  is  within  the  angle  of  collected  data.   However, 
the  polar  plots  for  the  higher  frequencies  show  increasingly 
asymmetrical  patterns  near  the  acoustic  axis  and  lack  of  a 
well-defined  lobe  pattern  expected  at  those  frequencies.   This 
non-uniformity  and  deviation  from  a  point  source  lobe  renders 
the  data  above  50  kHz  unreliable. 

At  the  25  cm.  range,  the  data  up  to  50  kHz  shows  close 
agreement  with  theory  for  the  low   angles  of  0°  and  10°. 
However,  as  the  angle  is  increased,  the  interference  from 
the  reflected  signal  causes  an  increasing  oscillation  of  the 
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data  about  the  theoretical  curve.   The  plots  for  the  30  cm. 
range  show  similar  results,  although  the  lower  angles  do 
not  compare  as  well  as  at  25  cm.   This  can  be  attributed 
to  the  concurrent  diffraction  signal  from  the  back  edge  of 
the  plate.   A  plywood  extension  was  added  to  the  plate  to 
see  if  the  data  at  30  cm.  could  be  improved.   The  extension, 
with  resulting  reduction  of  back  edge  diffraction,  improved 
the  data  slightly,  except  for  one  or  two  points,  which  may 
be  due  to  the  impedance  mismatch  between  the  plywood  and 
aluminum.   From  these  results,  the  geometry  was  restricted 
to  frequencies  no  greater  than  50  kHz,  grazing  angles  of  0° 
and  10°,  and  a  range  of  25  cm. 

C.   ROUGH  SURFACE 

The  backscatter  from  a  randomly  rough  surface  due  to  a 
transmitted  signal  at  normal  incidence  was  measured  in  order 
to  compare  to  experimental  and  theoretical  scattering  results 
reported  by  Clay  and  Medwin  [2].   Previous  studies  deal  with 
surfaces  with  small  slope  roughnesses  and  conclude  that  the 
backscattered  intensity  will  be  the  sum  of  a  coherent  com- 
ponent and  an  incoherent  component  with  the  coherent  dominant 
at  g  <<  1  and  incoherent  at  g  >>  1.   The  measured  statistic 
of  rms  height,  a,    should  be  comparable  to  the  value  calculated 
from  acoustical  normal  incidence  backscatter  in  the  range  of 
g  <  1  as  described  in  Section  III.B.   The  lead  shot  surface 
presented  a  backscattered  intensity  (not  normalized)  at  nor- 
mal incidence  that  was  almost  identical  in  shape  and  deviated 
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less  than  5  dB  from  the  backscattered  intensity  from  a  smooth 
surface  (Figure  25) .   The  fact  that  for  most  frequencies  the 
return  from  the  wedge  covered  with  lead  shot  yielded  a  slightly 
greater  signal  than  from  the  smooth  surface  may  have  been 
the  result  of  incoherent  contribution  of  backscatter  from 
multiple  large  slope  off-axis  scatterers.   The  similarity  of 
these  curves  and  the  small  dimensions  of  the  lead  shot  compared 
to  the  sound  wavelength  led  to  the  tentative  conclusion  that 
this  surface  appeared  nearly  "acoustically  smooth"  for  the 
frequency  range  of  interest.   The  term  "acoustically  smooth" 
is  used  in  the  classical  sense,  g  <<  1. 

The  experiment  now  shifted  to  data  collection  from  sur- 
faces coated  with  gravel  to  present  a  greater  degree  of 
roughness . 

Two  methods  were  used  to  derive  the  statistical  parameters 
of  the  surface.   In  the  acoustical  technique,  the  normal  inci- 
dence mean-squared  backscattered  pressure  (normalized  to 
squared  pressure  received  from  smooth  surface  at  the  same 
range)  was  plotted  as  in  Figure  26.   The  backscattered  pressure 
was  measured  at  four  regions  of  the  rough  surface,  to  obtain 
the  mean  squared  pressure.   Comparison  of  the  curve  with 
theory  (eqs.  27,  29)  was  used  to  extract  the  rms  height  and 
correlation  length.   Alternatively,  a  small  sample  plate  of 
gravel  roughness  was  constructed.   The  heights  in  a  5  cm.  square 
area  were  sampled  every  0.1  mm.  to  provide  a  2500  point  data 
base  for  calculating  the  same  statistical  parameters. 
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Theoretically,  for  g  <<  1,  the  log  of  the  relative  inten- 
sity is  proportional  to  the  roughness  parameter,  g.   The 
roughness  parameter  (g)  is  related  (see  eq.  22)  to  the  rms 
height  (a),  which  can  now  be  obtained.   The  latter  parameter 
and  the  asymptotic  value  that  the  curve  approaches  enable 

calculation  of  the  correlation  length  (L) .   The  backscatter 

2 
data  were  plotted  against  f   and  estimates  from  data  in  Figure 

26  were  calculated  as  follows, 


TBS  m         -g 


IBS 

10  log  ~£  =   -io  g  log  e 

X0 


g   =   k  a   (cos  9,  +  cos  Q -) 


2    2 

for  90°  backscatter,   g  =  4k  a 

I  2  2 

10  log  -^  =   -4.343  (16tt  f  )c2       c  =  345  m/s 
X0  c 

A  [10  log  =£2.] 

0  -7   2 
= =   slope   =  (-5.762  x  10   )a 

A[f^] 


a2   =   6.94  x  10"3  cm2 

-2 

a  ,  =  8.30  x  10    cm   =   0.83  mm 

acoustical 

a        =2.23  mm 
msrmt 
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and 


ASYMPTOTE 


AL2 

2~~ 2 
g>>l      2c   r 


r  (ASYMPTOTE)  2a  2r2  ,1/^2 


ASYMPTOTE  =  0.38 

-2 
a  =  8 . 3  x  10    cm 

r  =  25  cm 

2 
A  =  579.  cm   from  window  opening 


L, _   _..   ,  =  7.73  x  10  '  cm  =   0.77mm   using  a 

■  1C3  .L 


L..     ,  .   ,  =  2.0  mm  using  a„ 
Acoustical  ^   Msrmt 


L„     =4.0  mm. 
Msrmt 


The  statistical  measurements  performed  by  Jordan  [11]  on  the 
sample  gravel  plate  yielded  values  of  a    =   2.23  mm.  and  L  =  4.0 
mm. 

Although  the  shape  of  the  plotted  data  seems  promising, 
the  differences  in  the  two  statistical  parmaeters  are  evident 
by  the  plot.   The  slope  shown  leads  to  a  =  0.83  mm.  vice  surface 
measured  a  =  2.2  3  mm.   Similarly,  the  acoustically  determined 
correlation  length  is  0.77  mm.  vice  the  surface  measurements 
4.0  mm.   The  large  discrepancies  between  measured  and  acous- 
tically determined  statistical  parameters  can  be 
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explained  by  the  inappropriate  theoretical  description  of 
the  rough  surface  and  the  inability  to  use  a  large  number  of 

scattering  surfaces  for  the  experimental  determination  of 

2    2 

<P>    /P   •   The  experiment  used  gravel  that  presented  a  jagged 

rough  surface  with  steep  slopes  and  cusp-like  features, 
whereas  the  acoustic  scattering  theory  from  which  a  and  L  were 
calculated  assume  a  gently  changing  slope  for  the  rough 
surface. 

The  experiment  next  proceeded  to  identify  the  effects  of 
both  a  wedge  and  the  rough  surface  on  the  backscattered  sound. 
The  backscattered  signal  strength  (BSS  )  for  the  rough  surface 
alone  is  expressed  by  the  equation  for  incoherent  point 
scatters , 

IBS    R4 
BSSQ   =   10  log[-^2.  — £— ]  . 

0  R0   A 

This  definition  of  BSS  assumes  that  the  intensity  diverges 

-2  -2 

as  R   from  source  to  scatterer,  and  R   from  point  scatterer 

-4 
to  receiver  for  a  net  divergence  of  R   .   The  point  scattering 

assumption  was  verified  for  frequencies  above  50  kHz  at  nor- 
mal incidence  backscatter.   This  supports  the  assumption  of 
point  scattering  at  all  lower  frequencies,  and  at  grazing  as 
well  as  normal  incidence.   The  total  signal  returned  from  a 
rough  surface  and  wedge  is  the  sum  of  these  scattering  effects, 


13      1       4 

BSS   =   10  log[-^  -E, +  -M   ?   3  . 

0  rQ   X    10  RQ   A 
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The  experimental  data  is  a  total  intensity  ratio  from  the 
composite  surface,  where  the  proportion  due  to  the  individual 
effects  is  unknown.   The  geometry  of  0°  and  10°  at  25  cm. 
range  was  used  so  that  the  results  from  the  wedge  as  reported 
in  Section  B  could  be  separated  from  the  plots  with  confi- 
dence that  the  remaining  information  was  a  result  of  the 
rough  surface.   The  lack  of  any  specific  shape  of  the  data, 
Figure  27,  that  might  possibly  equate  to  a  simple  frequency 
dependence,  is  found  also  in  the  literature  [3,12];  as  the 
backscatter  due  to  sound  at  low  grazing  angles  is  determined 
by  complex  processes  of  multiple  scatter  and  shadowing  which 
have  not  been  adequately  described  by  any  theory.   The  ampli- 
tude of  this  near  grazing  backscatter  in  the  plots  is  signifi- 
cant in  that  it  completely  masks  any  diffraction  effect  of 
the  wedge,  which  is  5-15  dB  less  under  the  same  experimental 
conditions . 

This  absence  of  identifiable  wedge  backscatter  now  re- 
quired the  experiment  to  use  a  model  which  would  reduce  the 
rough  surface  backscatter  relative  to  wedge  diffraction. 
The  previous  conclusion  of  the  "smooth"  appearance  of  the  wedge 
with  the  lead  shot  covering  needed  to  be  reevaluated.   The 
slight  roughness  of  that  surface  at  low  grazing  angles  might 
produce  a  combined  backscatter  with  more  equitable  proportion 
between  the  magnitude  of  the  wedge  and  the  rough  surface 
return. 

The  low  grazing,  backscatterd  signal  from  the  lead  shot 
plate  was  sampled  over  different  time  apertures  to  see  the 
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varying  responses  of  rough  surface  backscatter  or  the  emer- 
gence of  the  diffracted  signal  from  the  crest  of  the  wedge. 
The  least  time  path  (LTP)  is  the  time  for  sound  to  travel 
directly  to  the  crest  of  the  wedge  and  return.   A  sampling 
window  prior  to  the  LTP  analyzes  scattering  from  the  rough 
surface  only;  whereas  a  window  that  includes  the  LTP  receives 
the  combined  scattering  effects  from  the  rough  surface  and 
the  wedge . 

An  additional  data  collection  run  obtained  a  greater 
frequency  resolution  by  changing  the  transmitted  frequency 
of  the  pulse  and  the  number  of  samples  taken  per  sampling 
window.   Figures  2  8  and  29  show  a  fairly  close  agreement 
between  the  slope  of  the  data  and  that  theoretical  slope 
for  the  wedge,  backscatter.   As  the  grazing  angle  becomes 
lower  (0°  situation,  Figure  29),  the  data  tends  to  diverge 
from  the  theoretical  line  as  the  magnitude  increases  from 
the  rough  surface  backscatter.   In  the  plot  of  increased 
frequency  resolution  (Figure  30) ,  the  sampled  aperture  of 
rough  surface  and  wedge  backscatter  (  •  )  follows  a  pattern 
of  linearly  decreasing  magnitude  with  increasing  frequency, 

yet  the  slope  is  much  greater  than  the  frequency  relationship 

1/2 
for  wedge  diffraction  alone,  which  is  theoretically   f 

[7].   In  Figures  32  and  33,  the  backscatter  from  the  rough 

surface  alone  presents  consistent  patterns;  although  once 

again,  no  specific  frequency  dependence  relationship  can  be 

modeled  from  the  data.   The  significant  feature  of  this 
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backscatter  effect  is  the  magnitude  of  the  rough  surface 
backscatter.   From  both  the  fact  that  magnitude  increases  as 
grazing  angle  decreases  and  the  relatively  greater  magnitude 
of  the  backscatter  as  compared  to  randomly  rough  surface 
backscatter  in  the  literature  [3,12],  the  results  suggest 
the  boundary  wave  effects  proposed  by  Tolstoy  [6] .   Tolstoy 
describes  a  boundary  wave  generated  from  a  rough  surface 
for  near  grazing  incidence  and  kh  <  1  (h  =  separation  of 
roughness  elements) .   The  magnitude  of  backscatter  for  the 
experiment  suggests  the  possibility  of  a  large  backscattered 
boundary  wave  for  surfaces  of  low  height  roughness,  a/A, 
but  large  slope  of  the  point  scattering  elements,  a/L,  for 
sound  at  low  grazing  angles. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

Experimental  analysis  of  impulse  wave  backscatter  has 
been  accomplished  in  the  laboratory  for  various  physical 
models . 

The  detected  return  from  the  shallow  change  of  slope  for 
a  plane  wedge  shows  close  agreement  with  the  theory  proposed 
by  Biot  and  Tolstoy. 

The  difference  between  the  actual  statistical  parameters 
of  the  rough  surface  and  those  estimated  from  the  spectrum 
of  scattered  energy  at  normal  incidence  is  attributed  to 
the  jagged,  cusp-like  appearance  of  the  surface;  that  is 
not  accounted  for  by  the  theoretical  description  of  the 
surface. 

The  inability  to  positively  separate  the  individual  back- 
scatter  identities  from  a  composite  surface  is  caused  by  the 
disproportion  in  magnitude  between  the  two  effects .   The 
total  magnitude  in  the  combined  spectrum  suggests  a  large 
backscattered  boundary  wave  for  low  grazing  angle  propaga- 
tion over  interacting  steep  sloped  rough  surface  elements. 

Additional  experimental  research  in  this  topic  could 
include  changes  or  improvements  in  the  areas  of  source/ 
receiver  selection,  physical  model  design,  and  theoretical 
description  of  rough  surfaces. 

The  ability  to  find  a  smaller  source/receiver  and  yet 
maintain  sufficient,  detectable  output  signal  and  receiver 
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sensitivity  could  extend  the  agreement  with  theory  above 
50  kHz.   Also,  the  design  of  an  electronic  switching  circuit 
to  allow  the  use  of  a  single  microphone  for  both  source 
and  receiver  would  expand  the  flexibility  in  geometry. 

A  rough  surface  wedge  with  a  greater  wedge  slope  would 
increase  the  backscatter  effects  from  the  wedge  so  that  they 
are  not  masked  by  the  return  from  the  rough  surface. 

Further  iteration  of  the  mathematical  description  of 
the  rough  surface  might  more  realistically  model  this  natural 
random  surface  to  the  experimental  data. 
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